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ABSTRACT 

The  generation  of  ambient  noise  by  physical  processes  dependent  on  shearing  of  the 
boundary  layer  under  sea  ice  is  investigated.  Special  attention  is  paid  to  the  identifica- 
tion of  individual  noise-generating  mechanisms  and  the  assessment  of  their  relative  im- 
portance. 

Recent  studies  of  Arctic  ambient  noise  are  reviewed  with  specific  reference  to  results 
showing  particularly  good  or  poor  correlation  between  ambient  noise  levels  and  ice 
movement  or  relative  current.  Potential  noise-generating  mechanisms  are  described  and 
categorized  according  to  their  small-scale  driving  forces  and  expected  noise  character- 
istics. More  detailed  quasi-objective  investigations  are  then  used  to  establish  the  relative 
importance  of  each  mechanism  as  a  contributor  to  the  overall  under-ice  noise  spectrum. 

Flow/Mechanical  mechanisms,  involving  ice  sheet  fracture  as  a  result  of  wind  and 
current-induced  bending  moments,  are  found  to  be  unlikely  contributors.  Conversely, 
processes  in  which  ice  fragments  in  current-driven  motion  under  the  ice  interact  to  cause 
bumping  and  grinding  noises,  appear  to  be  of  probable  importance.  Turbulent  pressure 
fluctuations  in  the  boundary  layer  under  sea  ice  are  shown  to  be  of  significance  at  low 
frequencies  on  a  local  scale.  The  role  of  resonant  cavities  in  the  under-surface  of  the  ice 
does  not  appear,  however,  to  be  an  important  one. 
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I.      INTRODUCTION. 

Recent  strategic  thinking  has  placed  great  emphasis  on  anti-submarine  warfare  in 
the  ice-covered  waters  of  the  Arctic  Ocean  (Young,  19S6).  This  emphasis  is  the  chief 
driving  force  behind  ever  increasing  efforts  to  fully  comprehend  and  model  underwater 
acoustics  in  this  region.  A  significant  proportion  of  this  effort  has  been  directed  towards 
the  study  of  ambient  noise  under  the  ice  and  in  the  marginal  ice  zone.  The  importance 
of  being  able  to  predict  ambient  noise  stems  from  its  variability  and  its  influence  on 
passive  sonar  performance. 

The  effectiveness  of  any  form  of  passive  sonar,  be  it  shipborne,  towed,  air  dropped 
or  mounted  on  the  sea  bed.  is  described  by  the  passive  sonar  equation.  From  Urick 
(1983) 

SL-  TL  =  XL-DI+DT 

where 

SL  is  the  source  level:  a  measure  of  the  sound  emitted  by  the  target 

TL  is  the  transmission  loss:  the  amount  of  attenuation  experienced  along  the  acoustic 
path  from  the  target  to  the  receiver 

NL  is  the  noise  level:    the  level  of  background  and  self  noise  obscuring  the  signal  at 
the  receiver 

DI  is  the  directivity  index:  a  measure  of  the  directivity  of  the  sonar  equipment 

DT  is  the  detection  threshold:  the  minimum  signal  level  that  can  be  detected  by  an 
operator  using  the  sonar  equipment. 

It  is  the  ability  to  accurately  measure  or  predict  the  magnitude  of  the  terms  in  this 
equation  which  permits  the  effective  determination  of  passive  sonar  performance.  A  firm 
knowledge  of  the  parameters  which  govern  the  passive  sonar  equation  therefore  allows 
an  operator  to  plan  both  offensive  and  defensive  operations  with  an  invaluable  degree 
of  foresight.  An  example  might  be  the  ability  to  forecast  a  time  period,  or  geographic 
location,  offering  poor  passive  sonar  conditions  in  which  a  submarine  might  safely  con- 
duct noisy  evolutions  while  avoiding  detection  by  hostile  forces.  Conversely,  a  unit  en- 
gaged in  offensive  operations  against  submarine  targets  under  the  ice  would  benefit  from 
an  ability  to  predict  a  good  acoustic  environment. 


The  magnitude  of  the  source  level  term  is  dependent  on  the  sound  energy  emitted 
by  the  object  of  the  passive  detection  effort.  It  is  therefore  governed  by  engineering, 
naval  architectural  and  tactical  considerations  not  directly  related  to  the  environment. 
Transmission  loss  depends  entirely  on  environmental  factors,  being  a  function  largely 
of  sea  water  and  sea  bed  characteristics.  Temperature  and  salinity  gradients,  both  hori- 
zontal and  vertical,  are  particularly  important.  Although  this  term  is  undoubtably  im- 
portant in  fixing  effective  sonar  ranges,  variability  over  moderate  distances  and  time 
periods  is  relatively  small  in  Arctic  waters.  This  means  that  good  climatological.  or  better 
still,  recently  measured,  profiles  of  temperature  and  salinity  are  usually  adequate  for  so- 
nar range  prediction  work  over  a  broad  area  and  for  considerable  time  periods. 

The  noise  level  term  is  generally  of  the  most  crucial  importance  in  the  overall  sonar 
performance  prediction  effort  in  Arctic  waters.  It  depends  to  a  large  degree  on  the  en- 
vironment and  is  known  to  be  highly  variable  both  spatially  and  temporally.  Noise  level 
incorporates  the  effects  of  self  noise,  created  by  the  hydrophone  itself  or  by  the  platform 
on  which  it  is  mounted,  and  ambient  noise  from  the  surrounding  water  mass.  Ambient 
noise  from  the  environment  is  a  variable  factor  depending  on  many  environmental  forc- 
ing functions.  Its  accurate  prediction  is  crucial  to  the  successful  assessment  of  passive 
sonar  performance. 

The  last  two  terms  of  the  equation,  directivity  index  and  detection  threshold,  are 
functions  of  equipment  design,  operating  technique  and  in  the  case  of  the  latter,  opera- 
tor effectiveness. 

Thus  it  can  be  seen  that  ambient  noise  prediction  is  an  essential  part  of  planning 
for  Arctic  operations.  An  additional  utility  can  also  be  found  in  weapons,  sensor,  and 
underwater  communications  systems  design,  where  a  grasp  of  typical  ambient  noise  lev- 
els is  important. 

Considerable  effort  has  been  directed  towards  the  modelling  and  prediction  of  am- 
bient noise  under  ice.  The  problem  has  been  tackled  in  a  variety  of  ways  including  em- 
pirical approaches,  attempts  to  predict  noise  levels  through  the  summation  of  emissions 
from  sources  and  by  consideration  of  energy  budgets  (Oard,  1987).  Pritchard  (1988) 
emphasised  the  need  to  investigate  and  understand  the  physics  behind  noise-generating 
processes  on  a  local  scale  before  modelling  them  and  summing  their  effects  at  a  partic- 
ular location.  The  overall  predicted  noise  would  then  be  the  sum  of  the  noise  intensities 
from  many  locations,  generated  by  a  variety  of  mechanisms  and  subject  to  transmission 
losses  depending  on  receiver  depth,  source  depth,  range  to  the  source  and  frequency.  The 
resulting   Arctic   ambient    noise   model   would    be    driven    by    outputs    from   existing 


meteorological  and  ice  dynamics  models.  A  categorization  scheme  for  noise-generating 
mechanisms  was  proposed,  by  Pritchard,  based  on  forcing  factors.  Mechanisms  impor- 
tant under  compact  ice  conditions  are  separated  into  three  areas 

1.  ridging  effects  which  are  driven  by  large  scale  ice  stress  and  deformation 

2.  micro-cracking  effects  which  are  driven  by  local  ice  stress 

3.  effects  dependent  on  the  shearing  of  the  boundary  layer  of  the  ocean  under  the  ice. 

It  is  the  purpose  of  this  study  to  identify  and  investigate  any  mechanisms  which  may 
fall  into  the  last  category.  Such  potential  noise-generating  mechanisms  are  catalogued 
as  exhaustively  as  possible  and  categorized  according  to  their  small  scale  driving  forces 
and  expected  noise  characteristics.  Investigations  are  then  conducted  with  a  view  to 
establishing  their  importance,  if  any,  in  the  overall  generation  of  ambient  noise  under 
ice. 

Many  of  the  physical  processes  involved  are  highly  complex,  and  some  are  poorly 
understood  or  the  subject  of  current  research.  For  this  reason  the  investigations  are  to 
some  degree  subjective.  An  attempt  has  been  made  to  model  the  physics  in  as  elementary 
a  fashion  as  possible  and  generally  to  draw  inferences  from  'order  of  magnitude'  results. 

Chapter  II  is  a  review  of  recent  Arctic  ambient  noise  research  which  considers  the 
effects  of  boundary  layer  shear.  The  purpose  of  this  brief  review  is  twofold.  Firstly  to 
explore  the  fundamental  concept  that  noise  is  caused  by  shearing  of  the  boundary  layer 
under  ice.  and  secondly,  to  identify  individual  mechanisms  which  have  been  suggested 
as  being  responsible  for  noise.  Evidence  is  sought  to  establish  whether  shearing  in  the 
oceanic  boundary  layer,  in  the  form  of  ice  movement  through  the  water  or  water  move- 
ment under  stationary  ice.  is  important  as  a  forcing  function  for  noise  generation. 

Chapter  III  outlines  potential  noise  generating  mechanisms  and  a  proposed  cat- 
egorization scheme.  The  individual  mechanisms  are  then  investigated  in  Chapter  IV  with 
conclusions  and  a  summary  presented  in  Chapter  V. 


II.     REVIEW  OF  CURRENT  LITERATURE. 

Previously  published  studies  of  under-ice  ambient  noise  are  fairly  numerous  and  in- 
vestigate a  broad  range  of  different  forcing  functions.  Forcing  by  thermal  effects  were 
investigated  by  Ganton  and  Milne  (1965)  and  Milne  (1972)  while  forcing  by  pressure 
ridging  and  ice  mechanics  was  studied  by  Pritchard  (1984)  and  by  Buck  and  Wilson 
(19S6).  Various  other  atmospheric  forcing  functions  have  also  been  studied  by  Payne 
(1964;.  Ganton  and  Milne  (1965)  and  Greene  and  Buck  (1978). 

Three  recent  papers,  Makris  and  Dyer  (1986),  Lewis  and  Denner  (1988)  and  Waddell 
and  Farmer  (1988).  consider  shearing  of  the  oceanic  boundary  layer  in  the  form  of  ice 
movement  or  flow  of  water  under  ice  as  a  forcing  function.  These  papers  are  therefore 
reviewed  briefly  with  the  objective  of  highlighting  results  which  tend  to  confirm  or  deny 
the  importance  of  boundary  layer  shearing  under  ice  as  a  forcing  function  for  noise 
generation. 

Makris  and  Dyer  (1986)  analysed  low  frequency  ambient  noise  measurements  and 
comprehensive  environmental  data  from  the  1982  FRAM  IV  expedition  to  the  central 
Arctic  (83°AT20°£).  The  octave  band  10  -  20  Hz  was  taken  as  being  representative  of 
noise  in  the  low  frequency  range  and  compared  with  temperature  and  various  applied 
stress  components.  The  underlying  hypothesis  was  that:- 

Noise  is  created  by  ice  fracture  mechanisms  proportional  to  the  state  of  stress  in  the 
ice  as  induced  by  environmental  loading. 

Cross-correlation  coefficients  between  the  low  frequency  noise  records  and  shear 
stresses  due  to  the  wind  and  current  were  found  to  be  particularly  high,  0.84.  in  both 
cases  for  9.9-day  records  and  0.71  and  0.74,  respectively,  for  23.7-day  records.  The  cor- 
responding coefficient  for  tensile  stress  due  to  cooling  was  0.15.  Composite  measure- 
ments, namely  ice  stress  on  the  ice  sheet's  vertical  section  and  the  stress  moment  acting 
about  the  ice  sheet's  central  horizontal  plane,  were  also  highly  correlated  with  ambient 
noise,  correlation  coefficients  being  0.81  and  0.87,  respectively,  for  9.9  day-records. 

The  general  conclusion  of  the  work  was  that  low  frequency  pack  ice  noise  correlates 
best  with  the  moment  due  to  opposing  wind  and  current  stresses  acting  on  the  ice  and 
worst  with  air  temperature. 


Lewis  and  Dcnner  (1988)  studied  the  correlation  between  noise  at  10,  32  and  1000 
Hz  and  a  number  of  parameters  describing  temperature  and  ice  kinematics.  The  data 
analysed  were  obtained  during  the  AIDJEX  project  in  the  Beaufort  Sea  during 
1975-1976  and  was  exceptionally  comprehensive.  High  quality  positional  data  were 
available  from  some  40  satellite-tracked  drifters  which  allowed  the  accurate  determi- 
nation of  ice  translational  speed  and  differential  kinematic  parameters.  Correlation  co- 
efficients were  calculated  for  ambient  noise  levels  at  the  above  frequencies  and  the 
following  parameters: 

U,  the  translation  rate  of  the  ice 

U2,  proportional  to  the  kinetic  energy  of  the  ice 

D,  the  divergence  of  the  ice 

|D| 

V,  the  vorticity  (rotation  rate)  of  the  ice 

I.V| 

N,  the  normal  deformation  rate  of  the  ice 

S.  the  shear  deformation  rate  of  the  ice 

F,  the  temperature  difference  between  the  air  and  the  ice 

and 
(X2  +  S:)1 2,  the  total  deformation  of  the  ice. 

The  results  of  these  calculations  are  reproduced  as  Figure  1.  The  correlation  ob- 
tained throughout  between  ambient  noise  and  L',  the  simple  translation  rate  of  the  ice 
measured  by  satellite  tracking,  is  quite  striking.  In  virtually  all  cases  the  correlation  was 
the  maximum  observed  for  any  of  the  parameters  above.  U  correlated  particularly  well 
with  32  Flz  noise  in  the  summer  (August  1975).  This  was  attributed  to  the  ice  "rushing 
through  the  water".  Another  observation  was  that,  in  general.  U  correlated  better  than 
U2,  except  at  1000  Hz.  This  suggests  that  pure  ice  motion  is  more  important  in  produc- 
ing noise  than  the  kinetic  energy  of  the  floes,  except  at  higher  frequencies  where  collision 
events  are  probably  more  significant. 
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Figure  1.        Linear  correlation  coefficients  between  ambient  noise  and  ice  kinematic 
parameters  (from  Lewis  and  Denner,  1988). 


Overall  this  work  seems  to  emphasize  the  importance  of  ice  movement  in  the  gen- 
eration of  ambient  noise.  Although  no  particular  mechanisms  are  mentioned,  the  role 
of  shearing  in  the  oceanic  boundary  layer  appears  to  be  of  probable  significance. 


Waddell  and  Farmer  (1988)  describe  the  noise  resulting  from  the  break-up  ofland 
fast  ice  in  Amundsen  Gulf  in  the  Canadian  Archipelago.  Although  the  study  was  in  fact 
motivated  by  the  possibility  that  ice  conditions  could  be  monitored  acoustically,  the 
conclusions  are  of  some  interest  here. 

Data  collected  over  a  127-day  period  from  April  to  August  1986  included  ambient 
noise  over  the  frequency  range  300  Hz  to  14.5  kHz,  current  measurements  and 
meteorological  observations.  Currents  in  the  area  are  strongly  influenced  by  the  tides. 
This  is  shown  in  Figure  2  which  also  depicts  time  series  of  longwave  infra-red,  ambient 
noise  and  air  temperature  for  a  10-day  period  in  June. 
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Figure  2.       Time  series  of  ambient   noise  and  various  environmental   parameters 
(Waddell  and  Farmer,  1988). 


Initially,  when  the  ice  was  consolidated  and  land  fast,  there  was  no  obvious  corre- 
lation between  the  current  meter  observations  and  ambient  noise.  A  more  marked  cor- 
respondence was  noted  between  noise  and  atmospheric  factors.  A  strong  northerly 
current,  developing  in  July,  however,  coincided  with  a  sharp  increase  in  noise  levels  as 
the  ice  near  the  recording  instruments  began  to  break  up  (see  Figure  3). 


60 


Figure  3. 


ambient  sound  record 


Time  series  showing  the  increase  in  ambient  noise  levels  with  increased 
current  and  ice  breakup   (Waddell  and  Farmer,  1988). 


These  results  seem  to  indicate  that  when  ice  is  landfast  and  static,  a  moderate  cur- 
rent under  it  does  not  create  discernable  noise  in  the  range  300Hz  to  14.5  kHz.  When 
the  ice  cover  is  in  motion,  however,  significant  correlations  can  be  found  between  cur- 
rent speeds  and  high  frequency  noise.  This  is  probably  be  due  to  collisions  and  cracking 
during  ice  breakup. 

As  a  general  summary,  the  studies  done  by  Makris  and  Dyer  (19S6)  support  the  idea 
that  low  frequency  (10-20  Hz)  noise  results  from  opposing  wind  and  current  stresses 
acting  on  the  ice.  The  results  of  Lewis  and  Denner  (19SS)  show  strong  correlation  be- 
tween noise  and  ice  translational  speed  which  is  proportional  to  the  relative  current  un- 
der the  ice.  Waddell  and  Farmer  showed  that  little  noise  results  from  relative  current 
under  stationary  ice. 


III.     POTENTIAL  NOISE  GENERATING  MECHANISMS. 

Ambient  noise  due  to  shearing  of  the  boundary  layer  under  sea  ice  could  be  caused 
by  a  number  of  different  mechanisms.  Little  or  nothing  has  been  written  about  the  de- 
tailed processes  involved  and  in  some  cases  the  physical  processes  are  sufficiently  com- 
plex to  prohibit  meaningful  analytical  investigation. 

A  number  of  potential  noise-generating  mechanisms  are  put  forward  here.  The  list 
is  not  intended  to  be  exhaustive  but  most  mechanisms  which  appear  to  have  any  possi- 
bility of  being  relevant  to  real  world  conditions  are  included.  Some  mechanisms  appear 
intuitively  more  important  than  others  while  some  may  seem  a  little  obscure.  All  are  are 
included  in  this  chapter  and  will  be  investigated  to  a  greater  or  lesser  degree  in  Chapter 
IV. 

A.     A  CATEGORIZATION  SCHEME. 

Two  broad  categories  of  potential  noise  generating  mechanisms  have  been  identi- 
fied. Firstly.  Flow  J  Mechanical  mechanisms  which  involve  mechanical  interaction  between 
the  sea  and  the  ice  cover  and  produce  sound  through  cracking  or  impactive  events. 
Secondly.  Flow  mechanisms  which  are  purely  a  function  of  fluid  flow  in  the  presence  of 
the  solid  ice  boundary. 

The  potential  mechanisms  are  therefore  divided  and  ordered  in  this  study  according 
to  these  criteria.  Figure  4  summarizes  the  mechanisms  that  have  been  considered  and 
their  relative  positions  in  this  scheme  of  organisation.  Figure  5  is  a  schematic  summary 
showing  the  partition  of  Flow; Mechanical  and  Flow  mechanisms. 
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Figure  4.      An  overview  of  potential  noise-producing  mechanisms. 
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Figure  5.      A  schematic  summary  showing  the  partition  of  Flow/Mechanical  and  Flow 
mechanisms. 
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B.     FLOW/MECHANICAL  MECHANISMS. 

1.  Cracking  Noise  Due  to  Wind/ Current- Induced  Moments. 

The  concept  of  opposing  wind  and  current  stresses  inducing  moments  in  the 
Arctic  ice  cover  was  originated  by  Makris  and  Dyer  (19S6).  Although  no  detailed  de- 
scription of  the  mechanism  itself  was  presented,  the  high  correlation  obtained  in  the 
FRAM  IV  data,  between  noise  and  the  stress  moment  acting  about  the  ice  sheet's  hori- 
zontal plane,  indicates  the  potential  importance  of  this  mechanism. 

In  an  attempt  to  investigate  the  detailed  mechanics  of  the  processes  involved, 
the  wind  and  current  stresses  are  taken  to  have  most  effect  on  features  with  some  ver- 
tical extent.  Pressure  ridges  or  ice  hummocks  present  obvious  points  where  winds  and 
currents  could  act  to  set  up  bending  moments  in  the  ice  sheet  which  might  possibly  lead 
to  ice  fracture  and  therefore  noise.  The  magnitudes  of  such  bending  moments  are  a 
function  not  only  of  the  relative  velocities  of  the  wind  and  under-ice  current,  but  also 
of  the  size  and  shape  of  the  ridge  and  of  the  thickness  and  general  nature  of  the  ice  sheet 
itself. 

Maximum  opposing  wind  and  current  velocities  are  most  likely  to  occur  when 
the  ice  is  in  a  low-stress  state  with  some  degree  of  wind  driven  motion.  Ice  in  a  static 
state  is  less  likely  to  be  subjected  to  strong  induced  moments  as  the  geostrophic  and  tidal 
currents  in  the  open  Arctic  tend  to  be  relatively  weak.  Wind  and  current  stresses  will 
also  be  greater  in  areas  of  highly  deformed  ice  cover  as  drag  there  is  more  pronounced. 
Fracture  as  a  result  of  induced  bending  moments  may  not  neccessarily  be  more  likely. 
however,  as  significant  deformation  is  often  coincident  with  greater  ice  thickness,  lower 
salinity  and  therefore  greater  strength. 

2.  Ice  Fragment  Bumping. 

The  importance  of  ice  collisions  in  the  production  of  ambient  noise  has  been 
recognized  in  several  works  (Diachok  and  Winokur  ,1974:  Shepard,  1979)  and  is  an  es- 
tablished mechanism  when  considered  on  a  large  scale,  i.e.,  floe-floe  interactions.  The 
presence  and  movement  of  significant  amounts  of  ice  debris  under  sheet  ice  is  also  doc- 
umented (Zubov,  1943;  Buck  and  Wilson,  1986).  A  potential  noise  generation  process 
therefore  exists  in  any  pack  ice  covered  region  where  deformation  events  create  quanti- 
ties of  loose  fragmented  material  which  can  be  moved  around  by  a  relative  current. 
Noise  from  this  source  is  likely  to  be  aurally  detectable  and  characterized  by  bumping 
or  grinding  sounds. 


13 


C.     FLOW  MECHANISMS. 

1.  Noise  From  Turbulence. 

Turbulence  is  a  recognized  source  of  low  frequency  noise  in  the  open  ocean 
(Wenz,  1962;  1972).  Although  the  processes  by  which  it  is  radiated  are  exremely  ineffi- 
cient and  its  far  field  effects  are  minimal,  the  direct  effect  on  hydrophones  within  tur- 
bulent flows  can  be  significant.  This  is  because  the  pressure  fluctuations  impinge  directly 
on  the  active  surfaces  of  hydrophones  which  react  by  generating  an  electrical  signal. 
Thus  turbulent  pressure  fluctuations  can  cause  a  type  of  self-noise  whose  effects  may 
be  significant  at  very  low  frequencies. 

In  the  Arctic  the  situation  is  complicated  by  the  presence  of  ice  cover.  The  ice 
acts  as  a  boundary  of  variable  roughness  past  which  water  flows.  In  addition,  the  ice 
sheet  itself  may  play  a  role  in  the  radiation  of  turbulence  noise  to  ranges  beyond  the 
actual  source  region.  Turbulent  self-noise  has  been  recognised  as  a  possible  contaminant 
of  Arctic  ambient  noise  data  by  Makris  and  Dyer  (19S6)  and  is  included  as  a  noise- 
generating  mechanism  in  its  own  right. 

Flow-induced  noise  and  flow-excited  noise  are  described  briefly  in  Chapter  IV 
as  potential  mechanisms  by  which  the  effects  of  turbulent  pressure  fluctuations  are  ra- 
diated to  the  far  field. 

2.  Resonant  Cavities. 

The  under-surface  of  floating  ice  in  Arctic  regions  is.  by  its  very  nature,  rough, 
deformed  and  pitted  with  cavities  caused  by  melting,  fracture  or  brine  rejection.  Such 
holes  or  cavities  present  a  potential  noise  source  when  water  is  in  motion  across  their 
openings.  Flow  may  be  the  result  of  ice  motion  or  currents  under  static  ice  and  could 
excite  resonant  responses  in  these  cavities. 

No  known  reference  has  been  made  to  the  possibility  that  resonant  cavities 
might  be  responsible  for  under-ice  noise  but  similar  above-water  features  are  certainly 
capable  of  producing  noise  in  strong  winds.  Noise  from  hull  cavity  resonance  is  also  well 
known  in  naval  architecture  (Blake,  1986). 

The  parameters  affecting  the  frequency  and  intensity  of  sound  from  resonant 
cavities  would  be  cavity  size,  shape  and  water  flow  speed. 
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IV.     AN  INVESTIGATION  OF  POTENTIAL  NOISE  MECHANISMS. 

A.     FLOW/MECHANICAL  MECHANISMS. 

1.     Wind/Current-Induced  Moments  at  Pressure  Ridges. 

The  action  of  a  current  relative  to  an  ice  field  on  submerged  pressure  ridge  keels 
and  the  action  of  the  wind  on  their  above  water  sails  gives  rise  to  a  potential 
flow  mechanical  noise-generation  mechanism.  The  moment  set  up  by  these  forces  may 
be  enough  to  cause  cracking  and  cracking  noise  in  the  ice  sheet  close  to  the  ridge  fea- 
tures. 

The  moment  Mk  due  to  water  pressure  on  the  ridge  keels  is  likely  to  be  rein- 
forced by  wind  action  on  corresponding  sail  formations  causing  a  moment  M,  (see 
Figure  6).  This  is  particularly  so  when  internal  stresses  within  the  local  ice  field  are  low 
and  motion  is  largely  due  to  wind-induced  free  drift.  The  applied  moment  M0  is  therefore 
a  function  of  the  wind  velocity  Vw  .  the  relative  current  velocity  Vc  and  the  pressure  ridge 
geometry.  Mg  is  the  sum  of  Mk  the  moment  caused  by  the  relative  current  on  the  keel 
and  M,  the  moment  on  the  sail  due  to  the  wind. 

This  functional  relationship  is  investigated  below  by  assuming  an  idealized 
pressure  ridge  model  after  Wright  et  al.  (197S)  and  Paquette  and  Bourke  (1988).  The 
geometry  of  this  model  is  illustrated  in  Figure  7. 

It  is  assumed  that 

1.  the  system  is  in  local  isostatic  equilibrium  when  V„  =  Vc  =  0 

2.  the  ridge  is  symmetrical  about  point  O  as  shown  in  Figure  7,  with  dimensions 
governed  by  local  isostatic  balance 

3.  the  density  of  the  ice  is  homogeneous  throughout  with  an  equal  percentage  of  void 
in  the  sail  and  keel  structures 

4.  the  surface  of  the  ice  sheet,  sail  and  keel  are  perfectly  smooth 

5.  the  wind  velocity  vw(z)  and  relative  current  velocity  vc(z)  vary  linearly  with  height 
and  depth  as  shown  in  Figure  8. 
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Figure  6.       The  action  of  the  wind  and  relative  current  on  a  pressure  ridge  feature. 
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Figure  7.       The  idealized  pressure  ridge  model. 
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Figure  8.       The  assumed  linear  wind/current  profile. 


a.     Ridge  Model  Geometry. 

From  the  assumption  of  isostatic  equilibrium 

y2hsivsPi  =  (hkivs+  Vihk(Wk-ws))(Pw-Pi) 


(i) 


where 
pw  is  the  water  density 
Pi  is  the  ice  density 


From  Figure  7 
and 

so 

2/u 

wk  =  — — +k;  4 

*       tan  aA.  ^  v  ' 

Substituting  (2)  and  (3)  into  (1)  and  rearranging 

=  0 


f   i    x 


K  , 

f      Pi     \ 

~ 2  h  " 

tan  a5 

K{ 

v   Pw-Pi   ) 

tan  c/.k 

h 

Solving  this  quadratic  in  [  — 


a  \        V  PW "  Pi  /  tan  a. 


^  /       Pi 


Pw  -  Pi 


(5) 


Now  taking  typical  values  for  p,  and  p„  as  920  kg  rrr3  and  1026  kg  m~3,  respectively,  and 
ak  =  30°,  ol:  =  20°  (from  Wright  et  al.,  1978) 

'h 
-p-^0.41  (6) 
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also  from  Fisure  7 


V  =  h  +  hk  +  —  tan  ak 


or 

tan  a.,. 

V-A»  +  *,-5E^  +  '*  <7> 

Since  the  ice  sheet  itself  can  also  be  assumed  isostatically  balanced 

hPi  =  'a(Pw  -  P/)  (8) 

and  since 
i  =  [s+  tk 

tk  =  t-£-*0.9t  (9) 

and 

t,-,H*p$-*0.L  (10) 

Substituting  (9)  into  (7) 

tan  y.k  pi 

hJ  =  hi,  +  /;, +  i  ——■ 

k         k       s  tan  o>s  Pw 

Further  substituting  numerical  values  for  cc,  a.k  and  hh  from  (6) 

hk'*Ahs  +  tk  (11) 

or 

hk'*\Jhk  +  tk  (12) 
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From  the  geometry  of  the  sail 
V  =  h5+ts 

and  from  (10) 

(pw  -  Pi) 


(13) 


h/  =  hs+  t 


Pv 


so 


hs'=*hs  +  0.\t 


(14) 


Now  that  functional  relationships  have  been  established  between  the  vari- 
ous ridge  dimensions,  an  equation  for  Ms,  the  moment  caused  by  wind  action  on  the  sail 
about  point  O  can  be  found  (see  Figure  9). 


Figure  9.       Wind  action  on  the  ridge  sail. 
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b.      The  Moment  Caused  by  Wind  Action  on  the  Sail. 

Using  concepts  from  Zubov  (1943)  and  the  theory  for  a  fluid  jet  acting  on 
a  flat  plate  (Hannah  &  Hillier.  1970)  the  wind  velocity  v„  at  height  z  above  sea  level  can 
be  resolved  into  two  components:  vp  the  component  parallel  to  the  surface  of  the  ridge 
sail  and  v„  normal  to  this  surface. 

vp  =  vwcosaJ  (15) 

vn  =  vw  sin  ccs  (16) 

At  the  windward  surface  of  the  sail  vn  causes  a  dynamic  force.  vp  has  no  effect  if  the  sur- 
face of  the  sail  is  perfectly  smooth.  The  resulting  force  P  is  equal  to  the  loss  of  momen- 
tum to  the  sail  per  second.    So 

—2 

P  =  rhvw  sin  as    kg  m  s 

where 
m  is  the  mass  of  fluid  (air)  incident  on  the  sail  at  height  z.  per  second. 

Then 

P  =  pavw  sin  0!v    kgms  (1 7) 

where 
pc  is  the  air  density. 

The  moment  m:{z)  about  point  O  is  given  by  the  product  of  the  length  of  the 
lever  arm  {  and  the  component  Pe,  of  P.  normal  to  i  (see  Figure  10). 
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Figure  10.      Moment  about  O,  the  point  of  symmetry,  due  to  the  wind. 


From  Figure  10 


/>,  =  />„  cos/? 


and 


p  =  it\l  -  6  +  a, 


so 


P{  =  P  sin(0  -  as) 
and  from  (17) 


2  . 


Pf  =  pavw  sin  as  sin(#  —  a.s) 


(18) 
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Now  sin(0  —  a,)  can  be  found  in  terms  of  h,',  /,  a,  and  z  by  considering  triangle  OAB,  (see 
Figure  11). 


Figure  11.      Triangle  OAB. 


From  Figure   1 1 


sin(#  —  Oj)  =  sin  6  cos  a5  —  cos  6  sin  as 


or 


sin(0  -  a5)  = 


x  cos  o.s  —  i  sin  as 


where 


x  = 


V  - ; 

tan  ar 


so 


sin(#  —  otj)  = 


/i/  cos  a5  —  i 
£  sin  otj 


(19) 
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Substituting  (19)  into  (18) 

PnK-(hr'  cos  aT  —  z) 
p,  =  fa  "l  I      1 (20) 

and 

ms(z)  =  SPt  =  paVl(hs'cos2<xs-z)  (21) 

If  the  wind  shear  is  linear,  as  outlined  in  the  assumptions  above 

>,(;)=  VJ:~'']  (22) 

"s 

so  substituting  (22)  into  (21) 

V2 
ms{z)  =  -^f-  (z  -  tf(h/  cos2*,  -  z)  (23) 

K 

The  total  moment  M,  on  the  sail  about  point  can  now  be  found  by  inte- 
grating with  respect  to  z  between  ts  and  h' 


o  V2    fV 

Ms  =  -^y^         (2  -  tf(hs'  cos2*,  -  z)<fc  (24) 


or 

K  =  ^TTT  Vh'\l  cos2*5  -  3)  +  V3/5(8  -  12  cos2*,)  +  hs'2t2(6  -  12  cos2*,)  -  Ahs't2s  cos2*,  +  ,54) 
1 2/C 
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From  (13) 

V  =  K  +  ts 

so 

V2 
Ms  =  -jf-  (h;(4  cos:a5  -3)  +  hsts(4  cos2a5  -  4)) 

or 

Ms  =  Pa\ili;Cs  (25) 

where  Cs  is  a  form  coefficient  dependent  on  the  geometry  of  the  ridge. 

(26) 

Ifoc^20°  and  i^OAr 

Q^(44-4^)xl0-3  (27) 

Figure  12  shows  Cs  plotted  against  hs  for  various  values  oft.  Note  that  C  ^  0.04  for  sail 
heights  greater  than  one  metre  and  is  virtually  independent  of  t.  If  y.s,  the  slope  of  the 
sail,  is  allowed  to  van",  an  interesting  effect  is  noticeable  (see  Figure  13).  Where 
af  >  28°  ,  Cs  reverses  sign,  i.e.,  the  result  is  a  moment  in  the  opposite  direction.  This 
could  be  relevant  when  considering  young,  relatively  unweathered  ridge  features  where 
slopes  may  be  greater  than  the  20°  assumed  above. 
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Figure  12.      The  form  coefficient  for  a  pressure  ridge  sail  as  a  function  of  sail  height 
and  ice  sheet  thickness. 


Figure  13.      The  form  coefficient  for  a  pressure  ridge  sail  as  a  function  of  sail  slope. 
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c.     The  Moment  Caused  by  Current  Action  on  the  Keel. 

The  moment  Mk  caused  by  the  relative  current  on  the  keel  can  be  found  in 
a  similar  way  to  that  found  for  the  sail.  The  relevant  geometry  is  outlined  in 
Figure  14. 


Figure  14.      Geometry  for  calculating  the  moment  caused  by  relative  current  on  a 
pressure  ridge  keel. 
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By  comparison  with  (21) 

mk  =  tPs  =  pwv?(hk'cos\-z)  (28) 

Linear  current  shear  is  assumed,  so 

*)-£^  (29) 

nk 

o   V2 

mk{z)  =  -^  (:  -  tkf(hk'  cos2a,  -  z)  (30) 

and 


v2  Ch 

Mk  - i^  (2  -  tk)\hk'  cos2ak  -  i)dz  (31) 

This  gives 

,-2 

Mk  =  ^tj-  hk(Ahk'  cos2**  -  Vik  -  4/,)  (32) 

From  (12)  h,'^\.lh-  +tk,  so 

yl 

Mk  =  -^y-  (/j2(6.8  cos2a&  -  3)  +  hktk{A  cos2a;,  -4))  (33) 

Ifas^  30°  and  /,  ^  0.9/ 

Mk—Z£f-  (2h2  -0.9V)  (34) 

or 

A4  =  Pvvlf/72Q  (35) 

where  Q  is  a  form  coefficient  for  the  keel  given  by 

163-75  j-  J  x  10"3  (36) 
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Figure  15  shows  Ck  plotted  against  hk  for  various  values  of/.  There  is  more  dependence 
on  i  than  in  the  case  of  Cs  and  it  should  be  noted  that  Ck ^ 0.15  for  hk>\  m.  Also 
comparing  Figures  12  and  15,  Ck  is  generally  about  four  times  bigger  than  C,  .  When 
ak  is  increased,  there  is  a  point  at  about  44°  where  Mk  reverses  sign  as  in  the  case  of  Ms 
at  greater  slopes  (see  Figure  16).  Again  this  may  be  relevant  when  considering  relatively 
new,  unweathered,  ridges. 
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Figure  15.      The  form  coefficient  for  a  pressure  ridge  keel  as  a  function  of  sail  height 
and  ice  sheet  thickness. 


Figure  16.      The  form  coefficient  for  a  pressure  ridge  keel  as  a  function  of  keel  slope. 
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The  total  moment  Ma  on  the  ridge  structure  about  point  O  is  given  by  the 
sum  of  the  moments  caused  by  wind  on  the  sail  and  relative  current  on  the  keel. 


Ma-MM  +  Mk 


(37) 


Figure  17  shows  Ma  with  its  constituent  components  M,  and  Mk,  as  given  by  (25)  and 
(35),  plotted  against  hs  with  t  =  1.5  m.  Typical  values  for  pfl,  VWt  pw  and  Vc  have  been  used, 
they  are  summarized  below: - 

pc=  1.3  kgm-3 

pw=  1026  kgm-3 

VK=  10  ms-1 

\\  =  0.2  ms-1 

These  wind  and  current  velocity  values  were  chosen  to  represent  the  upper  limits  of  ob- 
served measurements  from  the  literature  ,  i.e.,  a  worst  case  has  been  chosen  in  order  that 
the  resulting  moments  will  be  the  largest  that  could  normally  be  expected  to  occur  in 
nature. 


Mo. 

Mk 

Mt 


Vw  =  10m/ s    Vc  =  0.2m/s 


Figure   17.      Components  of  the  applied  moment  at  a  pressure  ridge  on  1.5m  thick  ice 
(kNm  per  m  ridge  length). 
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Figure  18  is  a  plot  of  Ma  against  t  and  ht  with  Kw=10ms-'  and  Vc  =  0.2  ms"1  .  The 
weak  dependence  on  ice  sheet  thickness  is  apparent  as  is  the  order  of  magnitude  of  M„ 
to  be  expected  for  these  conditions  ,  i.e.,  Ma  ~  1  kN  m  per  m  ridge  length. 
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Figure   18.      The  total  moment  applied  by  wind  and  current  at  a  pressure  ridge  as  a 
function  of  sail  height  and  ice  sheet  thickness  (kNm  per  m  ridge  length). 
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The  following  points  arise  from  the  preceding  calculation  of  the  moment 
due  to  wind  and  relative  current  action  on  a  pressure  ridge. 

1.  Mk  is  the  dominant  component  of  Ma. 

2.  t.  the  thickness  of  the  ice  sheet  has  little  effect  on  the  absolute  value  of  Mc. 

3.  For  angles  of  slope  o.s  and  a.k  taken  from  Wright  et  al.  (1978)  Ma  is  anticlockwise 
when  the  wind  blows  from  the  left  and  the  current  acts  from  the  right  (as  shown  in 
Figure  6).  If  steeper  slopes  are  assumed.  Mc  could  be  reversed.  Regardless  of  the 
sense  of  the  resulting  moment,  however,  the  absolute  value  of  Ma  at  the  ridge  will 
be  of  much  the  same  order. 

There  are  two  considerations  which  have  not  been  taken  into  account  when 
calculating  Mt  and  Mh  which  might  lead  to  increased  moments  at  the  ridge. 

1.  The  surfaces  of  both  the  sail  and  keel  have  been  assumed  smooth.  If  a  more  real- 
istic approach  were  to  be  taken  and  some  roughness  allowed,  drag  could  be  ex- 
pected parallel  to  the  sail  and  keel  surfaces. 

2.  Reduced  pressure  in  the  flow  downstream  of  both  the  sail  and  keel  would  also  have 
some  effect  on  M.  and  Mk.  This  would  contribute  in  a  clockwise  sense  for  a  wind 
blowing  from  the  left  and  current  from  the  right. 

Despite  these  possible  inaccuracies  and  any  others  which  might  arise  from 
the  assumptions  made,  equations  (25)  to  (28)  and  (35)  to  (37)  allow  a  reasonable  esti- 
mate of  Ma  to  be  made  for  sensible  values  of  ridge  height  and  ice  sheet  thickness. 
2.     The  Effect  Of  An  Applied  Bending  Moment  On  The  Ice  Sheet. 

In  order  for  the  moment  M„,  described  above,  to  cause  ambient  noise  under  the 
ice.  cracking  must  occur  in  the  ice  sheet.  The  horizontal  distribution  of  bending  moment 
M(x)  on  the  ice  sheet  must  therefore  be  investigated  in  order  to  estimate  the  maximum 
moment  that  could  be  expected  under  normal  circumstances.  Two  extreme  cases  are 
examined  to  determine  whether  the  maximum  expected  bending  moment  on  the  ice  sheet 
is  likely  to  be  sufficiently  large  to  cause  cracking.  The  first  case  examines  the  ridge  as 
a  relatively  new  feature  when  it  is  in  an  unconsolidated  state.  That  is,  it  is  assumed  that 
the  ridge  structure  consists  of  loose  rubble  and  ice  blocks  and  therefore  contributes  no 
strength  to  the  system.  The  second  case,  representing  an  older  pressure  ridge,  occurs 
when  the  ice  debris  comprising  the  sail  and  keel  has  frozen  into  a  consolidated  mass 
which  together  with  the  ice  sheet  near  the  ridge  can  be  considered  rigid.  Figure  19 
compares,  schematically,  the  two  conditions  outlined  above. 
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Figure  19.      Consolidated  and  unconsolidated  pressure  ridges. 
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a.      The  Unconsolidated  Pressure  Ridge. 

In  this  case  the  ridge  ice  sheet  system  is  assumed  to  behave  in  the  same 
way  as  two  adjoining,  loaded,  semi-infinite  beams  on  an  elastic  foundation.  The  uncon- 
solidated mass  of  rubble  and  ice  blocks  forming  the  sail  and  keel  interact  with  the  wind 
and  current  as  described  previously  to  cause  an  applied  moment  Ma  and  also  a  distrib- 
uted loading  on  the  ice  sheet.  They  do  not  however  contribute  in  any  way  to  the  strength 
of  the  sheet.  In  order  to  calculate  the  expected  bending  moments  on  the  ice  sheet  the 
section  to  the  right  of  the  ridge  axis  is  treated  as  a  semi-infinite  beam  with  static  loadings 
contributed  by  the  sail  and  keel  as  shown  in  Figure  20.  The  total  loading  can  be  thought 
of  as  being  due  to  three  triangularly  distributed  loadings.  That  is,  a  downward  load  due 
to  the  weight  of  the  ice  in  the  sail  (triangle  ABC  in  Figure  20)  and  an  upward  loading 
due  to  the  buoyancy  of  the  keel  (triangle  DEH  minus  triangle  FGH).  So  in  the  isostatic 
condition 

>'/(*)  =  }'aBC  -  XDEH  +  XFCH  (3S) 

where  y  is  the  downward  displacement  of  the  ice  sheet,  or 

y^x)  —y\dovn)  ~y\up)  09) 

Similarly 

Mi(.x)  =  MABC-MDEH  +  MrGH  (40) 

or 

Mi(.x)  =  M{down)-Mm  (41) 

M  is  the  bending  moment  on  the  ice  sheet  (positive  in  the  clockwise  direction). 
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Figure  20.      Static  loading  on  the  ice  sheet  by  an  unconsolidated  pressure  ridge. 
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Equations  for  the  displacement  and  bending  moment  on  a  semi-infinite 
beam  subject  to  triangularly  distributed  loading  have  been  developed  from  Hetenyi 
(1947)  in  the  appendix.  They  are:- 

y{x)  =  qQF{).Jtx)  0  <  x  <  if    (42) 

where  F(A,  f,  x)  =  (Q„_xl  +  AX{€  -  x)  +  BXxAu  -  DXx{Axe  +2BU)) 

y(x)  =  q0G(A,S,x)  x>€   (43) 

where  G{/.,  f.  x)  =  {Ci]e_xl  +  B)XAH  -  Dkx{A)t  +2BJ)  and 

M(x)  =  -q0H(X, f,x)  0<x   (44) 

where  H(/.,  {,  x)  =  {Ane_A  -  DixAu  -  Bix{Au  +  2B„)) 

{  is  the  base  length  of  the  triangular  loading 

and  q0  is  the  loading  at  the  free  end  of  the  beam,  i.e.,  at  x  =  0 

see  the  appendix  for  definitions  of  1//  the  characteristic  length,  k  the  foundation 
modulus  and  the  functions  Ax,  Bx.  C,  and  Dx. 

The  net  downward  loading  q0(dowr)  is  found  by  considering  triangles  ABC  and 
FGH  (see  Figure  20). 

%(down)  =  hsPiS  +  (/lk  ~  hk  -  tk)(pw  -  Pi)g  (45) 

and  using  (11)  and  (6) 

%{ftown)  =  hsS(Pi  +  l-6(pw  -  Pi)) 

Taking  p,  ^  920  kg  nr3  and  pw  ~  1026  kg  nr3 

?o(^vwo=10670/!s    kgnr's"2  (46) 

From  Figure  20 


3S 


For  the  upward  loading  qHup)  (due  to  triangle  DEH) 

%(uP)  =  (V  -  h)(Pw  ~  Pt)g  (4S) 

%(uP)  =  4hs(pw  -  pi)g 

or,  using  the  constant  values  given  above 

qom^4\55hs   kgnrV2  (49) 

From  Figure  20 

Substituting   for  q0  and  {  in  (42).  (43)   and  (44),   also  taking   k  =  (p„  -  p,)g  m  1039 
kg  nrV2 

}\down)(x)  =  °^42  F&>  * '(down),  x)  0<X<  { \dovn)    (51) 

y{down)(x)  =  0.942  G(A,S{down),x)  ^  *><W)    (52) 

*W)W  =  -4S5  //(/.,  e{down),  x)  0  <  .v    (53) 

>W)W  =  °- 144  F(;-  f(upy  x)  °^x<  em   (54) 

>W)W  =  °- J 44  GU >  r(upy  x)  x  ^  f(up)   (55) 

M{up){x)  =  -75  //(/,  e{up),  x)  0  <  x    (56) 
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Both  €{up)  and  €{dgwn)  are  related  to  hs  by  (47)  and  (50)  and  /  is  a  function  of  r.  so  from  (39) 
yfix)  =  0.942  F{down)(t,  hs,  x)  -  0. 144  F(up)(t.  hs,  x)  0<x<  t{down)    (57) 

y.{x)  =  0.942  G{down)(t,  hs,  x)  -  0. 144  F{up)(t,  hs,  x)  €{down)  <x<  f{up)    (58) 

yt{x)  -  0.942  C?(dkwfl)(rl  /V  x)  -  0. 144  G(w/J)(f,  ^,  x)  x  >  f{up)    (59) 

and  from  (41) 

Mix)  =  -485  H(down)U.  hs,x)  +  15  H(up)(t,  h5.x)  0<x   (60) 

Figures  21  and  22  show  the  displacements  of  ice  sheets  of  various  thick- 
nesses when  subjected  to  the  loading  of  unconsolidated  ridges  with  sail  heights  of  1.5m 
and  3m.  It  can  be  seen  that  the  concentrated  mass  of  the  sail  causes  downward  dis- 
placement at  the  ridge  axis.  The  more  distributed  buoyancy  force  due  to  the  keel  causes 
upward  displacement  away  from  the  ridge.  This  results  in  bending  moment  curves  as 
shown  by  Figures  23  and  24.  The  value  of  the  constant  /  was  calculated  using  a  typical 
value  for  E.  the  Young's  Modulus  for  ice.  £=3x  10"  Pa,  from  Mellor  (1986). 

The  absolute  value  of  the  bending  moment  due  to  the  isostatic  loading  on 
the  ice  sheet  reaches  a  maximum  value  .Vmax  some  10  to  20  m  from  the  ridge  axis.  It  is 
at  this  point,  with  or  without  the  additional  applied  moment  M„  due  to  the  wind  and 
relative  current,  that  cracking  is  most  likely  to  occur.  Figure  25  shows  .V/max  plotted 
against  hs  and  r. 

Having  established  estimates  for  Ma  the  applied  moment  and  .\/n,a,  the 
maximum  bending  moment  on  the  ice  sheet  due  to  isostatic  loading,  the  two  can  be 
compared  over  a  range  of  ridge  heights  and  ice  sheet  thicknesses.  From  Figure  18  and 
Figure  25  it  is  apparent  that  M0  is  of  the  order  0  to  1  kX  m  per  m  ridge  length  for  real- 
istic maximum  wind  and  current  speeds.  .Umax,  however,  is  of  the  order  20  to  2000  kN 
m  per  m  ridge  length,  i.e.,  .Wmax  is  roughly  two  orders  of  magnitude  greater  than  Ma  for 
a  given  ridge  size. 

Thus,  it  is  concluded  that  if  cracking  occurs  in  the  unconsolidated  case,  it 
will  be  due  almost  entirely  to  the  isostatic  loading  on  the  ice  sheet.  The  additional  ap- 
plied moment  M„  while  not  being  sufficient  to  cause  cracking  by  itself,  could  possibly 
act  as  a  trigger  mechanism  for  cracking  where  a  cracking  event  due  to  isostatic  loading 
was  likelv  to  occur  anvwav. 
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Figure  21.       Ice  sheet  displacement  due  to  isostatic  loading  of  an  unconsolidated 
pressure  ridge,  height  1.5m. 
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Figure  22.       Ice  sheet  displacement  due  to  isostatic  loading  of  an  unconsolidated 
pressure  ridge,  height  3.0m. 


41 


600 


500 


_400 


8300 


c20l 

x> 

c 

CO 


100 


V) 
JO 


o4^ 


/^•SStK^r...-., 


t  =  1.0  n 

t  =  1.5  ■ 

t  =  2.0  n 


10  20  30  40  50  60 

x  (metres) 


70  60  90  100 


Figure  23.      Absolute  bending  moment  on  the  ice  sheet  due  to  isostatic  loading  of  an 
unconsolidated  pressure  ridge,  height  1.5m. 
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Figure  24.      Absolute  bending  moment  on  the  ice  sheet  due  to  isostatic  loading  of  an 
unconsolidated  pressure  ridge,  height  3.0m. 
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In  order  to  establish  the  likelihood  of  cracking  due  to  isostatic  loading 
alone,  Mm3X  is  compared  to  the  critical  bending  moment  Mcnt  at  which  cracking  will  oc- 
cur. Mcrit  is  dependent  on  o,  the  tensile  strength  of  the  ice  sheet,  and  its  thickness  t.  A 
nominal  value  for  the  tensile  strength  of  the  ice  sheet  was  taken,  o,  =  0.4  x  106  Pa  from 
Mellor  (1986). 


M„lt  = 


<V 


(61) 


Figure  26  shows  the  ratio  MmJMe,„  plotted  against  h,  and  t.  From  this  graph  it  can  be 
seen  that  MmJMcril^\  where  t~h,.  When  MmJMcrit>\  cracking  due  to  isostatic  load- 
ing is  likely  to  occur.  It  should  be  noted  that  as  h,  becomes  much  greater  than  t  the  result 
is  a  greater  likelihood  of  cracking. 


Figure  25.      Maximum  bending  moment  on  the  ice  sheet  as  a  function  of  sail  height 
and  ice  sheet  thickness  (kN  in  per  m  ridge  length). 
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Figure  26.  The  ratio  between  maximum  bending  moment  and  critical  bending  mo- 
ment for  an  unconsolidated  ridge  under  isostatic  loading  as  a  function  of 
sail  height  and  ice  sheet  thickness. 
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b.       The  Consolidated  Pressure  Ridge. 

In  this  case  a  completely  rigid  system  is  envisaged.  It  is  assumed  that  the 
ridge  structure  has  undergone  repeated  freezing  events  and  is  no  longer  active  in  that 
ridge  building  processes  have  ceased.  The  material  comprising  the  sail  and  keel  has  fro- 
zen into  a  solid  mass  and  contributes  fully  to  the  strength  of  the  ice  sheet.  The  solid  ridge 
section  is  therefore  treated  as  a  rigid  beam  pivoting  about  point  O  (see  Figure  27)  while 
the  adjoining  ice  sheets  either  side  of  the  ridge  are  assumed  to  act  as  semi-infinite  beams 
on  elastic  foundations.  The  boundary  at  section  A-A,  between  the  rigid  portion  (the  sail, 
keel  and  an  unknown  length  of  ice  sheet)  and  the  rest  of  the  ice  sheet  (which  is  consid- 
ered flexible),  lies  a  distance  L  from  the  ridge  axis.  Equations  for  the  displacement  and 
slope  of  a  semi-infinite  beam  on  an  elastic  foundation  are  taken  from  Hetenyi  (1947)  and 
can  be  found  in  the  appendix. 

For  simplicity,  symmetry  about  the  ridge  axis  is  assumed  such  that 

1.  The  distances  LL  and  LR  equal. 

2.  yLandvR,  the  displacements  of  the  ice  sheet  at  distances  L  either  side  of  the  ridge 
axis,  are  equal  and  opposite. 

3.  ML  and  MR,  the  bending  moments  on  the  ice  sheet  at  distances  L  either  side  of  the 
ridge  axis,  are  equal  and  opposite. 

4.  QL  and  QR.  the  shear  forces  at  distances  L  either  side  of  the  ridge  axis,  are  equal  and 
opposite. 

Thus 

LL  =  —  LR  =  —  L 

yL  =  -yR=y 

ML  =  -MR  =  M 

and 

Ql  =  -Qr  =  Q 
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Figure  27.      Geometry  and  general  structure  of  the  consolidated  pressure  ridge  model. 
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Now  summing  the  moments  about  point  O 
LM0  =  Mb  +  ML  -  MR  -  LQL  -  LOR  -  Ma  =  0 

or  given  the  assumption  of  symmetry  above 

Mb  +  2M+2LQ  =  0  (62) 


Mb  is  the  moment  caused  by  the  immersion  of  ice  to  the  left  of  O  and  emmersion  of  ice 


to  the  rieht. 


Mb  =  I    x2k  tan  0  dx    =  2/3AL3  tan  0 

'-L 

or 

Mb^2l3kL36    for  small  6  (63) 

From  the  appendix,  the  displacement  of  a  semi-infinite  beam  subjected  to  a  shear  force 
O  and  moment  M  at  its  free  end  is  siven  bv 


2^_/}         2.U/2 
k        /x  k 


D/.X  I.        C?.x 


Sec  the  appendix  for  definitions  of  1//  the  characteristic  length,  k  the  foundation 
modulus  and  the  functions  Ax  BXt  Cx and  Dt  .  Since  the  section  A-A,  distance  L  from  the 
ridge  axis,  coincides  with  the  begining  of  the  flexible  ice  sheet  and  therefore  the  free  end 
of  the  semi-infinite  beam  model,  x  =  0  and  D>x  =  C>x  =  1.    Therefore 

20  A      ^u;2         2/ 
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The  slope  6  of  the  semi-infinite  beam  section  is  also  given  in  terms  of  Q  and  M  in  the 

appendix. 


AM/2 


A-    +  D 


k    •*** 

Again  at  the  free  end  x  —  0  so  Aix  =  Dix  =  1 
20  z1       4\/;3  ^;2 

e  =  --t~+/~~  =  -^W}--Q)  (65) 

Also,  from  Figure  27 

tan  0  =  y-  or  v  =  6L    for  small  0  (66) 

Now  subsituting  (63)  and  (65)  into  (62) 

4  L3,2(2M/.  -  g)  +2#n  +2LO  -  A/a  -  0  (67) 

and  from  (64).  (65)  and  (66) 

A//(l  +  2/.L) 

Q  =  — ; — ~, (68) 

*  (1+/L)  v    ; 

so  (67)  becomes 

M(  y  (/.L)3  +  4(/.L)2  +  4(/L)  +2)  -Mfl(l  +  (/.!))  =  0 

The  effect  of  the  applied  moment  Ma  on  M,  the  bending  moment  at  the  interface  region 
between  the  rigid  pressure  ridge  section  and  the  flexible  ice  sheet,  is  therefore  given  by 

-#--  <     ,"^;f" w 

iWfl     (y(;.Lr  +  4(;.L)2  +  4(;.L)  +  2) 

Figures  2S  and  29  show  the  ratio  MjMc  plotted  for  sail  heights  of  2  m  and  4  m  with 
various  ice  sheet  thicknesses.  /  has  been  calculated  using  a  Young's  Modulus  of 
£=  3  x  109  Pa  and  wind  and  current  speeds  of  10ms-1  and  0.2ms-1  have  again  been 
assumed  to  represent  a  worst  case. 


48 


0.5 


0.4 


0.3 


as: 


0.2 


0.1- 


2K 


\ 
\ 
\ 
\ 
\ 
s 

N 
\ 
\ 
S 
\ 


0.0 


t  =  2.0  ■ 

t  =  1.5a 

t  =  1.0  a 


■  i  i  »  i  i 1 

1O2O3O4O5O6O7D8090  10O 

L  (net res) 


Figure  28.      The  ratio  of  applied  moment  to  the  maximum  bending  moment  on  the  ice 
sheet  for  a  consolidated  pressure  ridge,  height  2m. 
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Figure  29.      The  ratio  of  applied  moment  to  the  maximum  bending  moment  on  the  ice 
sheet  for  a  consolidated  pressure  ridge,  height  4m. 
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Note  that  (69)  is  not  valid  for  L  <  —r-  since  the  keel  and/or  the  sail  struc- 
tures significantly  increase  the  rigidity  there.  It  is  assumed  that  cracking  is  unlikely 
within  the  ridge  structure  itself  due  to  its  increased  thickness  and  strength.  From 
Figure  29  the  maximum  bending  moment  on  the  ice  sheet  itself  appears  to  occur  where 
l        2     . 

The  maximum  bending  moment  Mmlx  is  plotted  for  various  values  of  sail 
height  and  ice  sheet  thickness  in  Figure  30.  Note  that  Mmtx  is  of  the  order  0  to  0.5  kN 
m  per  metre  of  ridge  length.  Comparing  this  with  Mcrl, ,  the  bending  moment  required  for 
cracking  to  occur  as  illustrated  by  Figure  31,  the  bending  moment  available  is  about 
three  orders  of  magnitude  smaller  than  that  required  for  cracking  to  occur  in  the  ice 
sheet. 

Thus,  it  is  concluded  that  cracking  noise  is  unlikely  to  occur  due  to 
wind/current  induced  moments  applied  to  a  consolidated  pressure  ridge. 


Figure  30.  The  maximum  bending  moment  on  the  ice  sheet  due  to  >vind/current- 
induced  moments  at  a  consolidated  pressure  ridge  (kN  in  per  m  ridge 
length). 
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Figure 


31.      The  critical  bending  moment  required  for  cracking  to  occur  in  ice  as  a 
function  of  thickness  (kN  m  per  m  ridge  length). 


3.     Conclusions  :  Cracking  Noise  Due  to  Wind/Current-Induced  Moments. 

From  the  results  of  a)  and  b)  above  it  is  apparent  that  the  generation  of  noise 
due  to  these  processes  is  unlikely  to  be  of  any  major  significance,  if  it  occurs  at  all  in 
nature.  In  the  case  of  the  young,  unconsolidated,  pressure  ridge  it  appears  that  the 
bending  moments  on  the  ice  sheet  due  to  the  isostatic  loading  of  the  sail  and  keel 
structures  far  outweigh  any  additional  bending  moments  caused  by  the  effects  of  wind 
and  current  on  the  ridge.  Any  cracking  which  occurs  is  likely  to  be  due  to  these  loadings. 
An  additional  contributor  to  local  ice  sheet  bending  might  also  be  large  scale  ice  stress. 
This  has  not  been  considered  here  as  it  would  probably  be  small  if  the  ice  were  in  wind 
driven  motion.  There  remains  the  remote  possibility  that  the  cracking  of  heavily  loaded 
ice  sheets  could  be  triggered  by  the  onset  of  wind/current-induced  effects.  It  seems  un- 
likely however  that  this  could  account  for  any  significant  correlation  between  wind, 
current  and  ambient  noise  in  field  measurements.  The  results  of  the  investigation  into 
the  more  mature  consolidated  pressure  ridge  reveals  even  less  probability  of  cracking  due 
to  wind/current  induced  bending  moments.  The  strength  of  the  ice  is  such  that  it  is  very 
unlikely  to  yield  to  the  relatively  weak  bending  moments  so  induced.  This  conclusion  is 
supported  by  the  fact  that  significant  noise  levels  have  never  been  observed  emanating 
from  old  inactive  pressure  ridge  features. 
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4.     Ice  Fragment  Bumping. 

This  mechanism  embraces  all  the  bumping  and  grinding  noises  which  would  be 
expected  to  occur  when  ice  fragments  move  about  under  an  ice  sheet.  The  probability 
that  this  mechanism  is  responsible  for  the  creation  of  significant  levels  of  ambient  noise 
is  reinforced  by  the  recorded  observation  of  audible  sounds  from  the  ice  surface.  Zubov 
(1943)  quotes  Weyprecht  as  stating  that 

Sometimes  shift  at  the  bottom  [  of  the  ridge  keel  ]  is  audible,  with  the  complete  re- 
pose of  the  ice  at  the  top.  This  occurs  probably  as  a  result  of  the  movement  of  water 
under  the  ice  field.  The  difference  of  the  movements  of  the  ice  field  and  of  the  water 
on  which  it  lies,  that  is,  the  current  of  the  water,  is  that  power  which  levels  the  lower 
depth  of  the  ice. 

Zubov  (1943)  also  quotes  Makarov  on  the  existence  and  movement  of  ice  blocks  under 

the  ice 

During  the  third  Winter,  a  fissure  was  formed  under  the  strip  under  the  Tram.  The 
fissure  opened,  lumps  of  significant  size  began  to  float  out  from  below.  "This 
shows",  says  Makarov,  "that  many  lower  lumps  constantly  travel.  The  current  of  the 
water  and  the  movement  of  the  ice  change  their  direction  so  that  the  migrating  ice 
block  stops  under  certain  conditions,  under  others  it  can  move  from  the  spot." 

It  is  neccessary  to  add  that  the  lumps  which  compose  the  under-ice  part  of  the 
hummocks  [  ridges  ]  are  not  only  transformed  but  are  also  gradually  destroyed. 

Ice  blocks  or  fragments  are  most  likely  to  be  the  products  of  active  ridging 
events.  In  the  previous  section  the  behaviour  of  young  unconsolidated  ridges  was  ex- 
amined quantitatively.  The  downward  depression  of  material  at  the  ridge  axis  was  shown 
in  Figures  21  and  22.  This  fractured  ice  appears  to  be  thrust  downward  and  trapped  in 
buoyant  equilibrium  under  the  ice  sheet  by  the  ridging  processes.  Buck  and  Wilson 
(19S6)  describe  the  forcing  of  ice  blocks  as  large  as  four  metres  thick  into  the  water 
column  in  their  study  of  an  active  pressure  ridge. 

It  is  reasonable  then  that  given  ice  sheet  motion  relative  to  the  water  body  be- 
low, the  loose  material  will  move  horizontally.  Indeed,  further  observations  by  Makarov. 
(Zubov,  1943)  suggest  that  this  is  the  process  by  which  pressure  ridge  keels  become 
vertically  truncated  and  more  widely  spread  in  the  horizontal  than  the  corresponding  sail 
formations. 
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Given  then  that  there  are  large  ice  fragments  in  motion  along  the  rough  under 
surface  of  an  ice  sheet,  collisions,  grinding  and  bumping  are  certain  to  occur.  Noise  from 
this  source  would  probably  be  similar  in  nature  to  that  observed  as  a  result  of  floe-floe 
interactions  in  the  marginal  ice  zone.  Predominant  frequencies  are  in  the  range  10  Hz  to 
1  KHz  according  to  Diachok  and  Winokur  (1974).  Certainly  the  impactive  nature  of  the 
processes  and  the  variety  of  fragment  sizes  is  likely  to  yield  noise  spread  over  a  wide 
frequency  range. 

The  intensity  and  perhaps  the  frequency  range  of  ambient  noise  caused  by  this 
mechanism  would  be  dependent  on 

1.  The  relative  current  under  the  ice. 

2.  The  proliferation  of  ice  debris  under  the  ice  sheet. 

3.  The  roughness  of  the  under-ice  surface. 

4.  The  thickness  and  density  of  the  ice  sheet. 

The  current  strength  is  of  prime  importance  and  is  directly  related  to  the  translational 
speed  of  the  ice  and  or  oceanic  or  tidal  currents.  The  presence  of  debris  and  ice 
roughness  are  both  a  function  of  the  deformation  history  of  the  ice  field  and  are  likely 
to  vary  by  region.  Item  4  has  been  included  as  these  factors  are  almost  certain  to  have 
a  bearing  on  the  nature  of  the  sound  resulting  from  impactive  events  and  its  propagation 
into  the  sea. 

Modelling  this  mechanism  is  likely  to  be  a  complex  task  if  approached  from  a 
purely  physical  standpoint.  A  more  profitable  approach  might  be  to  seek  empirical  re- 
lationships between  the  factors  outlined  above  and  measured  ambient  noise  levels.  This 
could  be  backed  up  by  some  simple  field  experiments  involving  perhaps  acoustic  moni- 
toring of  ice  blocks  pushed  through  bore  holes. 
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B.     FLOW  MECHANISMS. 

The  potential  noise  creation  mechanisms  falling  into  this  category  are  all  driven  ex- 
clusively by  the  interaction  of  the  sea  with  the  underside  of  the  ice  sheet  floating  at  the 
surface.  The  factors  important  in  their  operation  are  therefore  the  flow  velocity  of  the 
water  relative  to  the  ice,  the  characteristics  of  the  under-ice  surface  including  cavity 
shapes  and  sizes,  the  roughness,  and  to  some  extent  the  thickness,  of  the  ice  cover. 

It  seems  intuiuively  obvious  that  some  noise,  no  matter  how  small,  must  be 
produced  when  water  flows  under  rough  ice.  The  object  here  is  to  explore  the  mech- 
anisms which  may  be  responsible  for  this  flow  noise  and  assess  their  importance  in  the 
overall  generation  of  ambient  noise  in  an  ice  covered  environment. 

Because  of  the  complexity  of  the  mechanisms  outlined  in  this  section,  involving  as 
they  do.  turbulent  How  over  largely  random  surfaces,  the  prediction  of  the  noise  fields 
created  is  somewhat  crude.  Some  of  the  physical  processes  are  currently  being  researched 
and  many  of  the  concepts  are  still  only  poorly  understood  (Ross.  1987).  Where  possible, 
estimates  are  made  of  likely  sound  characteristics  in  terms  of  frequency  and  intensity. 
These  estimates  are  based  on  objective  considerations  but  must  be  assumed  accurate 
only  to  an  order  of  magnitude. 

A  typical  composite  Arctic  ambient  noise  spectrum  is  presented  in  Figure  32.  This 
serves  to  indicate  the  intensities  required  in  order  for  a  particular  mechanism  to  be  con- 
sidered important  as  a  contributor.  It  is  probable  that  (low  noise,  whether  or  not  it  is 
found  to  be  a  significant  component  of  the  total  noise  under  pack-ice.  would  be  entirely 
masked  by  noise  from  other  mechanisms  at  locations  near  the  ice  edge.  For  this  reason. 
the  measured  spectrum  chosen  to  represent  the  norm  is  from  a  central  arctic  pack-ice 
location.  It  is  reproduced  from  Makris  and  Dyer  (1986)  who  made  measurements  at 
S3\V20°£  in  April  19S2. 
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Figure  32.      A  typical  under-ice  noise  spectrum  (from  Makris  &  Dyer,  1986). 

1.     Noise  From  the  Turbulent  Boundary  Layer. 
a.      Boundary  Layer  Characteristics. 

The  ratio  of  kinetic  to  viscous  forces  in  a  flow  regime  is  described  by  the 
Revnolds  number  for  that  flow. 


Re  =  UQ  — 


(70) 


where 

U0  is  the  free  flow  velocity  beyond  the  influence  of  boundaries 

L  is  the  length  scale  associated  with  the  flow 

v  is  the  kinematic  viscosity  of  the  fluid  =  1.79  x  10-6nr2s  for  water  at  1  atmosphere 
and  0°C. 
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Whenever  the  Reynolds  number  exceeds  a  limiting  value  the  flow  becomes 
unstable  and  turbulence  results.  The  instability  is  caused  by  the  dominance  of  kinetic 
forces  in  the  flow  over  viscous  forces.  For  flow  under  ice,  where  relative  current  veloci- 
ties can  be  of  the  order  0.1  to  0.2  m  s  and  length  scales  measured  in  thousands  of  metres, 
Reynolds  numbers  can  easily  exceed  2000,  a  typical  limit  for  the  onset  of  turbulent  flow. 

In  the  presence  of  a  rigid  boundary,  which  in  the  first  approximation  can 
be  considered  flat  and  smooth,  a  turbulent  boundary  layer  is  formed.  The  thickness  5 
of  this  boundary  layer  is  a  function  of  the  Reynolds  number  for  the  flow  (White  ,1986) 
such  that 

Re 

where 

x  is  the  downstream  distance  from  the  begining  of  the  boundary  which  is  analagous 
to  the  length  scale  L  for  ice  sheet  cover. 

If  the  boundary  layer  under  the  ice  cover  is  fully  turbulent  it  can  be  con- 
sidered identical  to  the  mixed  layer  in  terms  of  constant  temperature  distribution.  Thus 
a  typical  value  of  30  m,  representative  of  arctic  mixed  layer  depths,  can  be  taken  for  S 
the  boundary  layer  depth.  It  is  accepted  that  processes  other  than  turbulence  in  the 
horizontal  flow  contribute  to  mixing  in  the  boundary  layer  under  ice  but  the  analogy  is 
presented  as  a  first  approximation.    From  (70)  and  (71) 

^-(^)i^r)";6  (72) 

=s  3000m  for  U0  =  0.2m  s-1  and  S  =  30m 

This  confirms  that  a  rough  length  scale  of  the  order  of  thousands  of  metres  with  a 
boundary  layer  depth  of  about  30  m  is  reasonable  for  large  smooth  ice  sheets  in  wind 
driven  motion. 
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b.     Noise  Mechanisms  in  the  Turbulent  Boundary  Layer. 

Vecchio  &  Wiley  (1973)  identified  two  mechanisms  by  which  noise  could 
be  radiated  from  a  turbulent  boundary  layer.  The  first  is  called  Flow-Induced  noise  and 
involves  acoustic  dipoles  induced  in  the  boundary  by  turbulent  pressure  fluctuations. 
The  second  is  called  Flow-Excited  noise  and  depends  on  the  elastic  behaviour  of  the 
boundary  which,  when  excited  by  turbulent  pressure  fluctuations,  radiates  energy  back 
into  the  fluid. 

The  direct  effect  of  the  pressure  fluctuations  themselves,  set  up  by  turbulent 
motion  in  the  boundary  layer  and  often  referred  to  as  Psuedo-Sound  or  Quasi-Sound. 
could  also  be  considered  as  a  noise  producing  mechanism  in  its  own  right.  Pseudo-sound 
is  a  non-radiating  phenomena,  it's  direct  effect  only  being  felt  when  pressure  fluctuations 
impinge  directly  on  the  active  face  of  a  hydrophone.  It  is  closely  related  to  self-noise  but 
is  considered  here  to  be  purely  a  product  of  the  turbulent  flow  and  not  dependent  on 
hydrophone  characteristics. 

Since  the  first  two  mechanisms  mentioned  above  depend  on  the  presence 
of  pseudo-sound  this  phenomena  will  be  examined  first,  both  as  a  driving  force  for  the 
other  mechanisms  and  as  a  source  of  noise  in  itself. 
2.     Pseudo-Sound. 

a.      The  Smooth  Boundary  Approximation. 

Skudrzyk  and  Haddlc  (1960)  describe  the  length  scales  of  turbulent  motion 
in  the  boundary  layer  as  varying  between  the  thickness  of  the  boundary  layer  and  the 
Kolmogorov  microscale  some  40  times  smaller.  Patches  of  turbulence  are  envisaged 
passing  over  a  hydrophone  or  presumably  any  given  point  on  the  boundary  at  a  rate 
Uc.  This  is  the  convection  velocity  which  is  proportional  to  U0  the  free  flow  velocity  such 
that  Uc  ^  O.8C0  .  The  characteristic  frequency^  of  the  resulting  pseudo-sound  pressure 
fluctuations  is  therefore  given  by 

/o=_^O.S—  (73) 

where 

d  is  the  separation  between  the  patches  of  turbulence  being  convected  past  a  given 
point. 
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If  d  is  approximately  equal  to  the  boundary  layer  depth  6 

Zo^O.8-^-  (74) 

Thus, /0^  0.007  Hz  for  U0  =  0.2  m  s"1  and  S=  30  m 

Skudrzyk  and  Haddle  (1960)  found  that  the  power  spectrum  of  pseudo- 
sound  pressure  per  unit  frequency  is  virtually  constant  for  frequencies  up  to  f0  and  then 
decreases  inversely  proportional  to  a  power  m,  at  frequencies  greater  than/;.  At  f=f0 
the  pressure  due  to  pseudo-sound  is  given  by 

pi  =  0.75  x  \0~5a2plulS*  4  {m  -  l/«0    Pa2  (75) 

where 
a  is  the  Kraichnan  constant 
L*0  is  the  free  stream  velocity 

<5*  is  the  displacement  thickness  of  the  boundary  layer  m  Sj5 
P„  is  the  water  density 

So  the  power  spectrum  for  pseudo-sound  is  given  by 

/>(/)  =  1.5  x  10-Vp2^  y  ('»  ~  Um)    Pa2Hz_1  /</0    (76) 

and 

P(f)  =  1.5  x  lO'Vp^d  -^  (m  -  1/m)  </0//f    Pa'Hz"1  f>f0    (77) 


a  is  related  to  the  Reynolds  number  for  the  flow  and  lies  between  about  0.7  and  6.  also 
Skudrzyk  &.  Haddle  (1960)  give  m  ^  3.   Therefore,  for  water  with  pw  =  1026  kg  nr3 

P(/)=1.6UoV<5    Pa2Hz"'  f<f0    (78) 

and 

P{f)  =\.6i-y d[f0!J)3    Pa2Hz_1  f>f0    (79) 
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Lotsch  (1971)  gives  an  alternate  formula  for  the  pseudo-sound  power  spec- 
trum in  terms  of  the  Strouhal  number  S 

POO  =  3.16  x  10~5 p2J*Ul{\  +  {nSf)-3'2    Pa2Hz_1  (80) 

where 

f5* 
5  =  ^Vt- 


'0 


or  takinc  <5*  =*  —  and  p„  =  1026  ks  m~3 
5 


(SI) 


Pressure  fluctuations  or  pseudo-sound  could  affect  a  hydrophone  located  in  the  turbu- 
lent boundary  layer  and  be  recorded  as  noise.  The  estimated  spectra,  as  given  by  (78) 
and  (79)  or  (SI),  when  compared  to  the  measured  spectrum  from  Makris  and  Dyer 
( 1 9 S 6 )  show  that  pseudo-sound  is  probably  insignificant  as  a  noise  source  if  the  under 
side  of  the  ice  is  smooth  (see  Figure  33).  Although  the  slope  of  the  roll-off  with  in- 
creased frequency  fits  well  with  the  measured  data,  the  overall  intensity  is  too  small 
compared  with  the  measured  spectrum  below  1  Hz. 
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Figure  33.        Comparison    of   the    calculated    pseudo-sound    specta    (assuming    no 
roughness)  with  a  measured  Arctic  noise  spectrum. 
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b.      The  Effect  of  Surface  Roughness. 

The  estimates  above  have  assumed  the  ice  sheet  to  be  a  smooth  flat  plate. 
This  is  a  poor  assumption  as  considerable  underside  roughness  is  known  to  exist 
(Wadhams.  19SS).  A  commonly  assumed  value  for  the  roughness  scale  under  ice  is  1 
m  although  higher  values  could  be  envisaged  in  some  circumstances. 

The  efTect  of  surface  roughness  in  boundary  layer  flow  is  to  pierce  the 
laminar  sub-layer  close  to  the  surface  of  the  boundary.  Vortices  are  shed  which  pass  into 
the  turbulent  region,  increasing  the  turbulent  motion  there  and  generally  shifting  noise 
levels  to  higher  frequencies.  These  effects  begin  to  be  important  when  the  Reynolds 
number  for  the  flow  at  the  roughness  scale  exceeds  a  value  of  about  5  (Skudrzyk  and 
Haddle,  1960),  i.e.,  when 

-^T— >5  (82) 


where 


U*  is  the  shear  velocity,  roughly  representative  of  the  velocity  in  the  immediate  vicin- 
ity of  the  laminar  sub-layer,  U*  a  O.O4c'0  (Skudrzyk  and  Haddle,  1960). 

and 

h  is  the  length  scale  for  the  roughness. 

For  h  =  1  m  and  [/0  =  0.2  m  sM 

-^f-  =*  5000  (83) 

Therefore  the  rough  underside  of  the  ice  can  be  expected  to  play  a  major  role  in  the 
creation  of  turbulence.  The  characteristic  frequency  for  pseudo-sound  created  by 
roughness  is  given  by 

O.Siq 

JQ{rough)  =         J  (S^) 

0.8L'0 

rather  than  the  — : —  given  by  (74)  for  the  flat  plate  approximation.  Thus  taking. 
L:0  =  0.2  m  s_1  and  h  =  1  m,  y^  2:  0.2  Hz  is  probably  a  more  reasonable  estimate  for 
pseudo-sound  under  rough  ice. 
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Equations  (78),  (79)  and  (SI)  can  now  be  rewritten  for  rough  boundary 
conditions.   The  Skudrzyk  and  Haddle  equations  (78)  and  (79)  become 

A/)=1.6y£/0Vc)    Pa'Hz"1  f<fQ{rouSh)    (85) 

and 

P{f)=\.6yUyd(f0{rough)IJ)2    Pa2Hz-]  f>fWough)    (86) 

The  Lotsch  equation  (81)  becomes 

/>(/)  =  6.1y5U3J\  +  (^)2j3!2    P^Hz"1  (87) 

where 

y  is  a  factor  to  take  into  account  the  increase  in  the  intensity  of  the  turbulent  pressure 
iluctuations  due  to  boundary  roushness. 
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The  next  problem  is  to  estimate  y.  If  the  noise  pressure  is  broadly  propor- 
tional to  the  drag  between  the  boundary  and  the  fluid  flowing  over  it,  as  suggested  by 
White  (1979)  y  will  be  approximately  equal  to  the  ratio  between  the  dimensionless  drag 
coefficient  in  the  rough  condition  and  that  in  the  smooth  condition.  White  (1979)  gives 
empirical  equations  for  the  dimensionless  drag  coefficient  CD  of  a  flat  plate  in  both 
conditions  and  also  presents  the  relationships  graphically  (reproduced  as  Figure  34). 
CD{smeBIh)  is  a  function  of  the  Reynolds  number  for  the  flow  such  that 

Q^)  =  ^  =  0.03l(4M-'/7  (88) 


A' 


c 


CDir0uSi:)  is  independent  of  Reynolds  number  in  the  fully  rough  condition,  depending  in- 
stead on  the  parameter  Lh  (see  Figure  34). 

trough)  =  (1-89  +  1.62  log  Lj/if512  (89) 

L  is  the  horizontal  length  scale  for  flow  over  the  plate 
h  is  the  roughness  scale 

An  estimate  of  the  increase  in  pseudo-sound  amplitude  for  a  rough  surface  will  therefore 
be  given  by  the  ratio 

(— 
m,L,h)=-E2^= L^_  (90) 

LD(smooih)       0.031(1.89+  1.62  log  UK)' 

Note  that  this  is  only  valid  for  fully  rough  conditions,  i.e.,  the  transition  y  -*  unity  for 
h  ->  0  is  not  described. 
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Figure  34.      Dimensionless  drag  coefficients  for  rough  and  smooth  fiat  plates  (from 
White,  1979). 
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Figure  35  shows  the  calculated  pseudo-sound  spectra  from  equations  (85) 
(86)  and  (87)  for  a  roughness  height  h  =  lm,  a  boundary  layer  depth  S  =  30  m  and  a 
free  flow  velocity  U0  =  0.2  in  s-1.  As  with  the  smooth  plate  approximation,  the  general 
shape  of  the  curves  fit  well  with  the  measured  spectrum  below  1  Hz.  The  characteristic 
frequency^  ,  however,  is  a  little  too  high  suggesting  that  the  effective  length  scale  may 
in  reality  lie  somewhere  between  3  and  h.  Despite  this  the  fit  looks  promising  given  the 
approximate  nature  of  the  theory  used,  and  the  fact  that  the  input  pararameters  are 
crude  estimations.  This  indicates  that  this  mechanism  is  probably  responsible  for  a  large 
proportion  of  the  noise  at  these  low  frequencies.  It  is  worth  noting  at  this  point  that 
Makris  and  Dyer  (19S6),  in  a  note  accompanying  their  composite  ambient  noise  spec- 
trum, hypothesise  that  the  noise  below  1  Hz  may  be  due  to  non-linear  surface  wave  noise 
or  pseudo-sound  in  the  oceanic  boundary  layer.  No  attempt  is  made  to  fit  a  curve  ex- 
actly to  the  measured  data,  as  none  of  the  free  variables,  i'Q  .  L  or  h  are  known  for  the 
data  set  used  to  produce  the  composite  spectrum. 

In  conclusion,  the  indications  are  that  pseudo-sound  may  well  be  an  im- 
portant source  of  noise  for  hydrophones  within  the  turbulent  boundary  layer  under 
rough  ice.  i.e..  positioned  no  more  than  a  few  tens  of  metres  below  the  ice.  The  fre- 
quencies involved  are.  however,  generally  below  1  Hz,  a  region  which  is  not.  at  present, 
of  general  importance. 
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Figure  35.      Comparison  of  the  calculated  pseudo-sound  specta  (roughness  effects  in- 
cluded) with  a  measured  Arctic  noise  spectrum. 
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3.  Flow-Induced  Noise. 

This  is  a  mechanism  by  which  the  effects  of  turbulent  pressure  fluctuations  are 
radiated  beyond  the  region  of  turbulence  to  the  far  field.  The  effect  is  due  to  two  phe- 
nomena. The  first  occurs  within  the  volume  of  the  turbulent  fluid  and  was  described  by 
Lighthill  (1954)  as  being  due  to  quadrapole  sources.  The  second  occurs  at  the  interface 
with  a  solid  boundary  and  is  described  as  being  dipole  in  nature  (Lotsch.  1971 ).  Acoustic 
power  due  to  the  first  of  these  effects  is  proportional  to  the  eighth  power  of  the  Vlach 
number  while  the  power  due  to  the  second  is  proportional  to  the  Mach  number  squared 
(Vecchio  and  Wiley,  1973).  In  the  under-ice  environment  Mach  numbers  are  very  small, 
typically  about  10~4 .  The  dipole  sources  therefore  dominate  but  the  total  radiated  power 
is  considerably  smaller  than  that  due  to  the  turbulent  pressure  fluctuations  themselves. 
Spherical  spreading  further  reduces  sound  intensities  due  to  this  radiative  mechanism. 
Flow-induced  noise  is  not.  therefore,  considered  to  be  an  important  contributor  to 
under-ice  noise. 

4.  Flow -Excited  noise. 

Vecchio  and  Wiley  (1973)  identified  this  mechanism  as  being  another  method 
by  which  noise  can  be  radiated  from  a  turbulent  boundary.  The  mechanism  depends  on 
the  presence  of  turbulent  pressure  fluctuations  along  a  flexible,  elastic  boundary. 
Vibrational  modes  are  excited  in  the  boundary  which  then  radiates  sound  back  into  the 
body  of  the  fluid.  Floating  ice  can  certainly  be  considered  as  a  flexible,  elastic  boundary; 
also  a  preceding  section  has  shown  that  low  frequency  pressure  fluctuations  are  created 
by  flow  over  its  rough  surface.  There  exists  a  possibility,  therefore,  that  this  effect  might 
be  responsible  for  noise  beyond  the  turbulent  boundary  layer  immediately  under  the  ice. 
The  physics  involved  are.  however,  extremely  complex.  For  this  reason  no  further  in- 
vestigation has  been  undertaken.  Detailed  analysis  of  this  process  is  warranted  and  may 
lead  to  the  conclusion  that  this  is  a  mechanism  of  some  significance. 
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5.     Resonant  Cavities. 

Water  flow  past  cavities  in  the  under  side  of  an  ice  sheet  is  a  potential  source 
of  noise  which  must  be  examined  in  order  to  establish  its  possible  importance  in  the 
overall  creation  of  ambient  noise.  Open  cavities  or  orifices  in  solid  boundaries  are  known 
to  resonate  under  certain  conditions.  The  onset  of  resonance  depends  on  the  fluid  flow 
velocity  past  the  cavity,  the  fluid  characteristics  and  the  shape  and  size  of  the  cavity  it- 
self. There  is  water  flow  past  cavities  in  the  under-surface  of  ice  sheets  or  the  submerged 
surfaces  of  pressure  ridge  keels  whenever  a  relative  current  exists.  The  object  of  this  in- 
vestigation is  to  determine  whether  this  flow  is  likely  to  be  sufficient  to  cause  noise  from 
such  cavities. 

a.      The  Helmholtz  Resonator  Model. 

The  Helmholtz  resonator  is  a  rigid  walled  cavity  whose  resonant  properties 
are  known.  Theory  pertaining  particularly  to  the  frequency  of  resonance  as  a  function 
of  cavity  dimensions  is  used  here  to  model  under-ice  cavities.  Figure  36  shows  cross- 
sections  of  the  classical  Helmholtz  resonator  and  an  idealized  simple  cylindrical  cavity 
such  as  might  exist  in  the  under  surface  of  an  ice  sheet  or  in  a  pressure  ridge  keel. 


Figure  36.      The  classical  Helmholtz  resonantor  and  a  simple  cavity. 
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F  is  the  volume  of  the  inner  cavity  which  provides  the  mass  element  of  the 
resonator.  S:  is  the  area  of  the  circular  opening  which  acts  as  a  simple  acoustical  source 
and  the  fluid  in  the  neck,  length  L,  provides  the  stiffness  or  spring  component.  From 
Kinslcr  et  al.  (1982)  the  frequency  of  resonance  is  given  by 

ft,»  =  c(i^)i;2  (9,) 

C  is  the  sound  speed  of  the  fluid 

L'  is  the  effective  length  of  the  neck  taking  radiation  -mass  loading  into  account.  For 
practical  purposes  L'  ^  L  +  Q.&A 

A  is  the  diameter  of  the  opening 

It  is  assumed  that  the  acoustic  wavelength  X  of  the  sound  resulting  from  resonance  is 
large  compared  with  the  length  scales  L.  Sl  2  and  K1 3.    For  the  idealized  ice  cavity 

L  =  B--j-  (92) 

S  =  J4r\  (93) 


and 

'•  =  T<f)3  CO 

Therefore 


w0=  C 


15 


1QAB-3A' 
Now  taking  C  ^  1500  m  s-1  for  sea  water 


co0^5S10(10/*Z?-3.42r'/2 


and 


fQ^925(\0AB-3AT      Hz  (95) 
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The  constraint  on  X  relative  to  the  cavity  dimensions  dictates  that  A  ~  D  for  (95)  to  be 
valid.    Figure  37  shows  resonant  frequencies  for  such  simple  cavities. 


Figure  37.      Resonant  frequencies  (Hz)  for  simple  cavities. 


b.     Excitation  of  Resonance. 

The  next  concern  is  whether  the  modest  flow  velocities  present  under  ice  are 
sufficient  to  excite  resonance  in  the  cavities  there.  From  Blake  (1984)  the  relationship 
between  the  free-flow  velocity  across  the  entrance  to  a  cavity  and/  the  frequency  of  the 
pressure  fluctuations  caused  by  the  flow  is  given  by 


feT-  =  -J^{n-\jA)     n=  1,2,3,.... 

U0  is  the  free  flow  fluid  velocity 

Cr  is  the  hydrodynamic  phase  velocity  across  the  opening 

A  is  the  diameter  of  the  opening 

C 

For  turbulent  boundary  layer  excitations  ——-  ^  0.33  ,  so 

0.33^0 
fe*  r-^in-l/A)     «=  1,2,3 


(96) 


(97) 
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Figure  38  shows  excitation  frequencies  fe  for  cavities  up  to  one  metre  in  diameter.  A 
nominal  value  of  U0  =  0.2  m  s-1  has  been  taken  to  represent  a  typical  maximum  relative 
current  velocity  under  ice. 
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Figure  38.        Excitation  frequencies  for  simple  cavities  (Hz)  in  a  flow  velocity  of  0.2 
ni/s. 

By  comparing  Figures  37  and  38  it  can  be  seen  that  the  frequencies  at  which 
under-ice  cavities  could  be  excited  by  a  typical  flow  velocity  are  far  lower  than  their 
resonant  frequencies. 

Thus,  it  is  concluded  that,  within  the  limits  of  the  approximations  made, 
resonant  cavities  are  unlikely  to  be  important  in  the  generation  of  ambient  noise  under 
ice. 
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V.     SUMMARY  AND  CONCLUSIONS 

The  importance  of  being  able  to  accurately  predict  ambient  noise  levels  in  Arctic 
waters  was  discussed  in  Chapter  I.  In  order  to  attain  that  ability  research  must  be 
undertaken  into  even-  aspect  of  the  Arctic  ambient  noise  problem.  The  reviews  of 
Chapter  II  and  the  statements  made  by  Pritchard  (1988)  illustrate  the  great  number  of 
different  subsidiary  problems  which  must  be  understood  and  overcome  before  a  reason- 
able solution  can  be  reached.  The  large  number  of  possible  mechanisms  involved  in  noise 
generation  and  their  extreme  complexity  make  this  a  very  daunting  task.  In  comparison 
to  the  problem  of  modelling  open  ocean  noise,  the  Arctic  problem  is  complicated  by  the 
different  interactive  processes  between  air.  ice  and  sea.  In  addition,  the  inaccessability 
of  the  Polar  regions  leads  to  a  dearth  of  meaningful  data.  Progress  is.  however,  being 
made.  Wenz  (1962)  commented  on  the  inadequacy  of  data  from  open  ocean  sites  at  that 
time,  but  was  able  to  describe  and  catalogue  the  noise-generating  processes  to  good  ef- 
fect. A  similar  situation  exists  now  with  regard  to  noise  produced  in  ice-covered  waters. 

Pritchard  (1988)  categorized  Arctic  noise  mechanisms  into  three  types  by  forcing 
function  (as  described  in  Chapter  I).  This  work  addresses  the  mechanisms  constituting 
the  third  category,  i.e.,  those  mechanisms  dependent  on  the  shearing  of  the  boundary 
layer  of  the  ocean  under  the  ice.  A  scheme  is  proposed  in  Chapter  III  by  which  this 
category  could  be  further  sub-divided  into  Flow  and  Flow  Mechanical  mechanisms.  The 
following  is  a  synopsis  of  the  results  of  investigations  into  this  sub-set  of  potential  noise 
mechanisms. 

A.     FLOW/MECHANICAL  MECHANISMS. 

1.     Cracking  Noise  Due  to  Wind/ Current- Induced  Moments. 

From  statistical  work  done  with  observed  data,  Makris  and  Dyer  (1986)  found 
that  ambient  noise  correlated  well  with  the  stress  moment  acting  about  an  ice  sheet's 
central  horizontal  plane.  In  this  study,  bending  moments  induced  in  the  ice  sheet  were 
investigated  by  considering  the  effects  of  wind  and  current  on  an  idealized  pressure  ridge 
feature.  Ridges  were  envisaged  as  being  in  one  of  two  conditions,  the  young  unconsol- 
idated ridge  and  the  more  mature  consolidated  feature.  In  the  case  of  the  first  it  was 
found  that  the  moments  caused  by  wind  and  current  are  insignificant  when  compared 
to  the  bending  moments  caused  by  isostatic  loading.  The  relationship  between  the 
bending  moments  present  due  to  isostatic  loading  and  the  critical  bending  moment  re- 
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quired  for  cracking  in  the  ice  sheet  was  identified  and  unstable  combinations  of  ridge 
height  and  ice  thickness  were  found.  There  appeared  to  be  a  remote  possibility  that  the 
wind  and  current-induced  bending  moments  might  act  as  a  precipitative  mechanism  for 
cracking  actually  resulting  from  isostatic  loading.  It  is  felt,  however,  that  this  is  not  a 
likely  explaination  for  the  correlation  observed  by  Makris  and  Dyer  (1986).  The  study 
of  the  consolidated  ridge  revealed  even  less  likelihood  of  wind  current-induced  cracking 
through  induced  bending  moments.  The  assumed  rigidity  of  the  system  resulted  in 
bending  moments  that  were  far  below  the  critical  value  required  for  initial  cracking. 
2.     Ice  Fragment  Bumping 

Beyond  some  observations  that  ice  debris  does  seem  to  move  about  in  current- 
driven  motion  causing  audible  sounds,  very  little  appears  to  have  been  written  on  this 
subject.  The  concept  is  attractive  in  that  similar  ice-ice  interactions  are  known  to  cause 
a  great  deal  of  noise  at  the  ice  edge.  Given  sufficient  amounts  of  freely  moving  fragments 
and  a  moderate  relative  current,  a  very  noisy  result  could  be  envisioned.  Although  the 
detailed  physics  of  this  processes  were  not  investigated,  there  are  indications  that  the 
resulting  sound  would  be  spread  over  a  wide  frequency  range,  possibly  varying  from  tens 
of  Hz  to  the  low  kHz.  This  is  a  promising  mechanism  warranting  further  investigation. 
The  problem  of  modelling  this  effect  in  practical  terms,  however,  does  present  many 
difficulties.  The  most  accute  would  probably  be  the  determination  of  the  density  of  free 
ice  fragments  in  any  given  region. 

B.     FLOW  MECHANISMS 

1.     Noise  From  Turbulence 

This  is  the  process  by  which  noise  is  produced  as  ice  moves  in  water  or  when 
water  flows  under  static  ice.  It  is  probably  the  process  envisioned  by  Lewis  and  Denner 
(19S8)  when  they  wrote  of  noise  caused  by  ice  "rushing  through  the  water."  Investi- 
gations of  the  turbulent  pressure  fluctuations  which  could  be  expected  for  typical  ice 
roughnesses  and  flow  regimes  revealed  a  predominance  of  very  low  frequency  activity. 
When  compared  to  a  measured  Arctic  spectrum,  the  predicted  turbulent  pressure  or 
pseudo-sound  spectra  have  the  correct  characteristics  to  explain  low  frequency  noise. 
Thus  it  is  suspected  that  pseudo-sound  in  the  boundary  layer  is  significant  at  frequencies 
below  1  Hz  but  only  for  receivers  actually  within  the  region  of  turbulence. 

Two  mechanisms,  Flow-Induced  and  Flow-Excited  noise,  by  which  the  effects 
of  turbulent  pressure  fluctuations  might  be  radiated  were  mentioned  briefly.  The  first 
appears  to  be  too  inefficient  to  radiate  significant  sound  energy.  The  second  may,  how- 
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ever,  be  of  some  importance.    Both  phenomena  need  further,  more  detailed,  consider- 
ation. 

2.     Resonant  Cavities 

The  idea  that  cavities  in  the  under  surface  of  the  ice  might  resonate  due  to  water 
flow  comes  from  a  consideration  of  similar  features  in  air.  A  comparison  with  the  clas- 
sical Helmholtz  resonator  concept  revealed  the  likely  frequencies  of  resonance  for  vari- 
ous cavity  dimensions.  By  establishing  the  frequencies  that  typical  flow  velocities  might 
excite,  however,  it  was  found  that  such  flows  are  almost  certainly  too  slow  to  cause 
resonance.  Thus  this  mechanism  is  considered  to  be  an  unlikely  candidate  for  noise 
production. 
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APPENDIX       BEAMS  ON  ELASTIC  FOUNDATIONS. 

This  appendix  describes  the  theoretical  behaviour  of  infinite  and  semi-infinite  beams 
on  elastic  foundations  as  given  by  Hetenyi  (1946).  This  theory  is  used  to  model  floating 
ice  because  the  ice  is  flexible  to  a  degree  and  since  buoyancy  and  gravity  make  the  sea 
act  like  an  elastic  foundation. 

A.      INFINITE  BEAMS. 

The  spatial  distribution  of  displacement.  y(x),  slope  6(x),  bending  moment  M{x)  and 
shear  force  O(x)  for  an  infinite  beam,  on  an  elastic  foundation,  subjected  to  a  point 
loading  P  and  moment  M  {see  Figure  39)  are  given  by:- 

y(x)  =  ^-A?x  +  -^f-B).x  (,41) 

d(x)  =  -^f-B/x  +  -^-C,x  (A2) 

M(x)  =  ^-C^  +  -fD,x  (,13) 

<?U)  =  --f  D/x--^-A/x  (AA) 

where 


A- is  the  foundation  modulus    =  (pw  —  pt)g  for  floating  ice  {AS) 

1//.  is  the  characteristic  length  of  the  system  =  (  —f—  V4    m-1  {A6) 

£is  the  Youngs  modulus  for  the  material  of  the  beam  per  unit  width 

I  is  the  moment  of  inertia  for  the  beam    =-r^-  (Al) 

i  is  the  thickness  of  the  beam. 
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-Ajc 


ALx  =  e      (  cos  Ax  +  sin  Ax) 


B?x  —  e  /JC  sin  Ax 


C?JC  =  e    '  (  cos  Ax  —  sin  Ax) 


D>JL  =  e    '"*  COS  ^.JC 


(AS) 

(A9) 

(AIO) 

(All) 


Figure  39.      The  infinite  beam  with  a  point  load  and  bending  moment. 

B.     INFINITE  BEAMS  WITH  TRIANGULARLY  DISTRIBUTED  LOADINGS. 

For  a  triangularly  distributed  loading  (see  Figure  40) 


<7o 


y{x)  =  ~aJo7  (c'-i'-*i "  c**  ~  1}{D>*  +  A){e  ~  x)) 


<Jo 


y{x)  =  !u7  (c^-*i  ~  c*  ~ 2UD^ 


M(x)  = 


% 


3       (A?U_x]-A^-2AJB;jc) 


% 


Q{X)  =  ~~4&  {B'U-xl  +  Bjjc  ~  ^^ 


where 
<70  is  the  load  at  x  =  0 
{  is  the  base  length  of  triangular  loading 


0^jt<<?  (All) 

x>{  (/H3) 

jc>0  (,414) 

;r;>0  (,415) 
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Figure  40.      Infinite  beam  with  triangularly  distributed  loading. 

C.     SEMI-INFINITE  BEAMS  WITH  TRIANGULARLY  DISTRIBUTED 
LOADINGS. 

If  the  infinite  beam  equations,  (Al)  to  (A4),  are  to  be  used  to  describe  a  semi-infinite 
beam  with  a  triangular  loading,  end  conditioning  forces  must  be  applied  in  order  to 
create  the  effect  of  a  free  end.  Thus,  if  a  free  end  is  to  be  created  at  x  =  0,  a  force  P'  and 
moment  M'  must  be  applied  to  exactly  cancel  Q0  and  M0.  The  result  will  then  be 
A/(0)  =  0  and  Q(0)  =  0  ,  the  conditions  for  a  free  end.    From  (A14)  and  (A15)  at  x  =  0 


A/0  = 


8;.v 


0H6) 


<2o  = 


% 


4/V 


{Bu-Xf) 


{All) 
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Substituting  (A16)  and  (A  17)  into  (A3)  and  (A4),  with  x  =  0 

M(0)  =  ^-  +  ^  +  -^7(^-l)  =  0  (,418) 

Q{0)  =  ^--M^---^-{BM-)I)  =  0  {A\9) 


Solving  (A  18)  and  (A  19)  for  F  and  A/' 

/>'  =  -^-  (  -^  +  1  -  2B}£  +  2).€)  {A20) 

2/ c 

M'  =  -  -^j  (  -A?/  +  1  -  Bki  +  ).{)  (.421) 

The  distribution  of  displacement  on  the  semi-infmite  beam  with  a  triangular  loading 
can  now  be  found  by  superimposing  the  effect  of  the  load  from  (A  12),  together  with  the 
effects  of  P'  and  M'  at  x  =  0  as  described  by  (Al). 

y{x)  =  ~xtl  (C*I'-*I  ~  C**  "  U€D**  +  AX{{  ~X))  +  ~2^  A**  +  "f"  B*   °  -  -T  <  {  {A22) 

>'W  =  "^7  (Q|^|  -  Q*  -  V./Dj^  +  -^f  ^  +  ■&-  B/x  x>f  {All) 

Substituting  for  P'  and  M'  and  simplifying 

&*>  =  4^7" (C'-"-*'  +  4/(^  ~'v)  +  ^*  "  D>AA>£  +  2B}e))  °  " v <  { {A2^ 

y{x)  =  ~mt  (c^-*i  +  ******  ~  D>*{A>*  +  2B"])  x~£  {A25) 

cPv 

And  since  M(x)  =  —  El—fr  (from  Hetenvi) 

dx- 

M{x)  =  -  TT7  {A>  i  '-*  i  ~  D***u  ~  B*AA*  +  2B^  x  -  ° ("* 26) 
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D.     SEMI-INFINITE  BEAMS  WITH  POINT  LOADINGS  AND  BENDING 
MOMENTS. 

Developments  similar  to  those  describing  the  behaviour  of  a  semi-infinite  beam  with 
a  triangularly  distributed  loading  are  used  in  the  case  of  a  point  loading  and  bending 
moment  at  the  free  end  (see  Figure  41). 


Figure  41.      Semi-infinite  beam  with  a  point  loading  and  bending  moment. 

Equations  for  the  displacement,  slope,  bending  moment  and  shear  force  distibutions 
are  given  below. 


,  v      IP).  n        2MT  r 
>M  =  -£—  D?JC —  C;jc 


lor  41/ j3 


(All) 
(A2S) 


MM—fBto  +  MAte 


Q{x)  =  -PCix-2M).BJix 


[A29) 
(A  30) 
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